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Suspensions are of wide interest and form the 
basis for many smart fluids [T143]- For most sus- 
pensions, the viscosity decreases with increasing 
shear rate, i.e. they shear thin. Few are reported 
to do the opposite, i.e. shear thicken, despite the 
longstanding expectation that shear thickening 
is a generic type of suspension behavior [9|. 
Here we resolve this apparent contradiction. We 
demonstrate that shear thickening can be masked 
by a yield stress and can be recovered when the 
yield stress is decreased below a threshold. We 
show the generality of this argument and quan- 
tify the threshold in rheology experiments where 
we control yield stresses arising from a variety of 
sources, such as attractions from particle surface 
interactions, induced dipoles from applied elec- 
tric and magnetic fields, as well as confinement of 
hard particles at high packing fractions. These 
findings open up possibilities for the design of 
smart suspensions that combine shear thickening 
with electro- or magnetorheological response. 

Shear thickening is presumed to be due to general 
mechanisms such as hydrodynamics 0, [13] or dilation 
[lll - [T3j . and thus all suspensions are expected to exhibit 
shear thickening under the right conditions 0. So far, 
however, the exact conditions have not been determined. 
One condition is apparently set by attractive particle in- 
teractions. It has long been known that attractions, ob- 
served as fiocculation in suspensions, can prevent shear 
thickening. This has been shown by modifying the chem- 
istry, for example by adding flocculating agents to ob- 
serve the transition from shear thickening to thinning (for 
a review, see@). In other cases, crossing the gel transi- 
tion was shown to eliminate shear thickening [TJ, Il5j . 
A key problem, therefore, is to understand how inter- 
particle attractions interfere with shear thickening. We 
demonstrate here that a simple and direct criterion for 
the existence of an observable shear thickening regime 
in dense, non-Brownian suspensions can be developed by 
comparing the yield stress produced by attractions with 
the inherent shear thickening stresses. We then gener- 
alize this condition to show how a yield stress from any 
source modifies the shear thickening phase diagram. 

Our experiments used an Anton Paar rheometer to 
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FIG. 1: Revealing shear thickening by adding surfac- 
tant to hydrophobic glass spheres in water. Soda lime 
glass spheres of mean diameter 90 /J,m with a hydrophobic 
silane coating were suspended at a packing fraction (f> = 0.52. 
a, A: viscosity curve without surfactant. The divergence of 
the curve is characteristic of a yield stress. •: viscosity curve 
of the same system at the same cf> with added surfactant. The 
shear thickening regime is the region of positive slope in the 
curves of viscosity n versus applied stress r. Shear thinning 
is characterized by a negative slope and Newtonian fluids, 
such as water, exhibit constant n. b, Images show clustering 
due to interparticle attractions (top) and no clustering when 
surfactant is added (bottom). Scale bar is 200 /im. All im- 
ages (including subsequent figures) were taken at rest under 
an optical microscope in a dilute quasi two-dimensional layer. 
In this dilute case, attractions can be observed by the high 
number of particle contacts in the form of clusters or chains. 



measure the shear stress r and the shear rate 7 of a 
wide range of different suspensions. The viscosity is 
defined as rj = r/7 in the steady state. Our focus is 
on non-Brownian, dense suspensions that show strongly 
packing-fraction-dependent, reversible shear thickening, 
often called 'discontinuous,' because of the steep stress- 
shear rate relationship. To understand the significance 
of interparticle attractions, we first consider the particle- 
liquid surface tension. Figure II Al shows the striking 
change in behavior produced by adding a small amount 
of surfactant to a water suspension of glass spheres with 
a hydrophobic coating. In the aqueous environment the 
coating leads to network-like particle clusters (Fig. II Ab . 
top) which minimize exposed surface area and thus po- 
tential energy. As a consequence, to pull particles apart 
requires overcoming a significant stress threshold. In 
Fig. II Al this translates into a region where, for applied 
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stresses smaller than this yield stress, the suspension docs 
not flow and the viscosity effectively diverges. The result 
is shear thinning behavior (A). Added surfactant elim- 
inates the clustering with its associated yield stress and 
reveals a region of underlying shear thickening (•) be- 
low the range of the previous yield stress. This suggests 
the yield stress due to attractions is responsible for hid- 
ing shear thickening if it overwhelms the shear thickening 
stress range. 

For a more detailed exploration than afforded by chem- 
ical means of the role played by the magnitude of the 
yield stress in modifying the shear thickening regime, we 
require in-situ, tunable control over the strength of the 
attractions. This can be achieved by applied electric and 
magnetic fields that polarize particles of a given dielec- 
tric or magnetic susceptibility and also have the practical 
advantage of reversibility. The result is a field-dependent 
attraction between neighboring particles and thus a con- 
tinuously tunable yield stress. We used dielectric glass 
spheres in mineral oil for electrorheology and magnetite- 
filled polyethylene glycol (PEG) rods suspended in PEG 
for magnetorheology. Figure II Dl shows the evolution of 
the yield stress and shear thickening regime with both 
types of field. For small fields, the viscosity curve is 
seemingly unaffected. A main result from these data 
is that increasing the field strength, and the concomi- 
tant yield stress, pushes the onset of shear thickening 
to higher stress values. At intermediate field values the 
curves rejoin the zero-field shear thickening behavior af- 
ter exhibiting a viscosity minimum. A yield stress thus 
simply results in a smaller range of applied stresses over 
which shear thickening is observable. Only when the 
yield stress becomes large enough that it encroaches on 
the upper limit of the shear thickening range is the effect 
fully eliminated. Qualitatively this behavior is neither 
dependent on the suspension nor the source of the yield 
stress, as seen from the similarity between panels a and 
c in Fig. II Dl The fact that the same conclusion can ap- 
ply to Fig. II Al is especially remarkable considering that 
the induced dipoles are directional but the chemical at- 
tractions are anisotropic. Given the different microstruc- 
tures, this indicates that it is the stress scale resulting 
from attractions that determines whether shear thicken- 
ing is observable or not. 

The experiments discussed so far concerned yield 
stresses produced by particle attractions. Similar behav- 
ior carries over to suspensions without attractive interac- 
tions in which a yield stress arises due to confinement at 
large packing fractions [l?} . Data are shown in Fig. II Al 
for several different packing fractions of cornstarch in wa- 
ter. It is seen again that the shear thickening range de- 
creases as the yield stress increases and eventually disap- 
pears when this yield stress approaches the upper stress 
limit of the shear thickening regime. 

The interplay between yield stress and shear thicken- 
ing emerging from the data in Figs. II Dl and II Al can be 
summarized in a set of non-equilibrium phase diagrams 
(Fig. II Dp showing the shear thickening, shear thinning, 



and jammed (defined here as a non-flowing state be- 
low the yield stress) regimes. Despite the differences in 
sources of a yield stress, there are important similarities. 
The stress thresholds bounding this regime (horizontal 
black lines) are nearly independent of <fi when the yield 
stress is small enough [HI [13, EH ■ As the yield stress 
increases, the lower threshold moves upward and even- 
tually approaches the upper boundary, at which point 
shear thickening ceases. For intermediate values of B, E, 
or 0, both jamming and shear thickening can be found 
at different stress values [IH HH HBJ . 

Since the boundaries of the shear thickening region are 
determined by local extrema of viscosity curves, they can 
be calculated given the relation between stress and shear 
rate in the lower shear thinning and shear thickening 
regimes. Note that the yield stress value is below the 
shear thickening phase boundary, leaving a shear thin- 
ning regime between the jammed and shear thickening 
regions. To quantify the effect of the yield stress on the 
shear thickening phase boundary, we therefore must ac- 
count for this additional regime. To model these contri- 
butions, we use the Herschcl-Bulklcy form, with a fixed 
exponent of 1/2 commonly used to describe shear thin- 
ning behavior [l(J 

T HB (i) =r y + a 1 j 1 / 2 . (1) 

Here the first term r y denotes the yield stress and a\ 
parameterizes the additional stress that is operative in 
the shear thinning regime. We refer to tub &s the shear 
thinning stress. Earlier work [lj| Ojl, [2(| suggests that 
contributions to the overall shear stress can be linearly 
separated as 

r(i) = T ffB (7) + a 2 7 1/£ , (2) 

where the second term represents the shear thickening 
stress parameterized by a 2 and an exponent e that ap- 
proaches zero in the limit where the stress/shear-rate re- 
lation becomes discontinuous at high packing fractions 
[l7| . Over the whole range explored in our experiments 
Eq. II Fl fits the data well, as demonstrated by the example 
in Fig. [5j 

The lower boundary of the shear thickening region oc- 
curs at the stress r m and shear rate j m corresponding 
to the local viscosity minimum. Differentiating ?/ = r/7 
and eliminating a 2 via Eq. II Fl evaluated at r m gives 

T m = T HB {j m ) + 2 ^ £ _ ^ [THBiim) + T y ] . (3) 

Eq. II Fl is in a form that directly shows how the shrink- 
age of the shear thickening regime depends on the shear 
thinning terms. The model parameters e, r y and a\ are 
obtained by fitting the data to Eq. II Fl for each value of 
B, E, and <j). The values of and e are found to be 
independent of B and E (see Suppl. Mat.), which can be 
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FIG. 2: Using magnetic and electric fields to tune the interplay between shear thickening and the yield stress. 

a, Viscosity curves for a suspension of ferromagnetic rods (254 x 32 x 25 ^m) for different values of applied magnetic field B. 
Magnetite-doped (30% by weight), PEG rods made by the PRINT process [16j were suspended in PEG at a packing fraction 
4> — 0.20. The shear thickening region is seen to shrink and eventually becomes eliminated as it is encroached on by the 
increasing yield stress, b, Microscope images show the rods for B — 30 mT (top) and B — (bottom), c, Viscosity curves for 
a suspension of dielectric spheres for different values of applied electric field E. Soda-lime glass spheres of diameter 90 /im were 
suspended in 58 mPa-s mineral oil at a cj> = 0.56. d, The microscope image shows the spheres for E — 60 V/mm. At E = 0, 
the image is similar to the bottom panel of Fig. II Ab . In both panels b and d the fields were applied vertically, in the direction 
of the shear gradient in a parallel plate rheometer. The scale bars are each 200 /im. 
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FIG. 3: Elimination of shear thickening by increasing 
packing fraction, a, Viscosity curves for cornstarch in wa- 
ter at different packing fractions <j>. The cornstarch particles 
had an average diameter of 14 /im. The water was density- 
matched to 1.59 g/mL by dissolving CsCl in it. A solvent 
trap was used to avoid evaporation and a Couette geometry 
was used to ensure the sample remained confined, b, The 
microscope image shows that particles do not cluster without 
confinement, also confirmed by optical tweezer measurements. 
Scale bar is 20 (im. 



seen from the overlay of shear thickening curves at higher 
stresses in Fig. IIP! indicating that the shear thickening 
stress is independent of attractions. This is in contrast 
to weaker, 'continuous' shear thickening due to hydrody- 
namics where attractions were found to affect the shear 
thickening stress [l5| . Equation II Fl is evaluated at the 
measured j m , T y , and a\ for each B, E, and <j), and a 
fixed value of e for each panel, as shown by the solid 
red symbols in Fig. II Dl This describes the lower phase 



boundary very well, typically within 12%. We note that 
Eq. II Fl along with an equation for j m (obtained from 
Eqs. II DIlIF! see Suppl. Mat.) can also be used to predict 
the phase boundary with attractions without measuring 
7 m for all field values, assuming only that the shear thick- 
ening stress is independent of the shear thinning mecha- 
nism. 

The agreement of Eq. II Fl with the measured phase 
boundaries in Fig. II D I demonstrates that the lower shear 
thickening phase boundary is set by the mechanism that 
produce shear thinning. Because of the second term in 
Eq. II Dl this is true whether or not the shear thinning in- 
cludes a yield stress. It is also independent of whether the 
shear thickening stress term is affected by the parameter 
that controls shear thinning if the shear thickening and 
shear thinning terms add linearly as seen in Fig. II Db 
and for 'continuous' shear thickening d, GJ]. The fact 
that this model reproduces the measured phase bound- 
ary confirms that, for 'discontinuous' shear thickening, 
the effect of attractions is to increase the shear thinning 
stress which hides shear thickening, rather than to af- 
fect the shear thickening stress directly. When B, E, or 
<p become large enough that shear thickening is not ob- 
served, r m (open red symbols in Fig. lID[) becomes higher 
than the shear thickening stress range. In the limit of 
e = 0, Eq. II Fl reduces to r m = thb(7 to ), so the stress at 
the phase boundary becomes equal to the shear thinning 
stress. Thus, the shear thickening regime starts to shrink 
when the shear thinning stress exceeds the stress at the 
onset of shear thickening, and it is eliminated when the 
shear thinning stress exceeds the stress at the viscosity 
maximum. Regardless of the particulars of the mathc- 
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FIG. 4: Non-equilibrium phase diagrams delineating observable shear thickening regions in terms of the associ- 
ated stress range. Stress range as a function of applied magnetic field B (a), applied electric field E (b), and packing fraction 
4> (c). The boundaries of the shear thickening regime are set by the local minima (T) and maxima (a) of the viscosity curves in 
Figs. II Dl and ll Al (•): the yield stress t v , below which suspensions are jammed. (•): the predicted onset of shear thickening r m 
evaluated from Eq. II Fl at the measured j m , demonstrating that the boundary is determined by the total shear thinning stress 
thb, regardless of the source of the yield stress. For panels a and b the values of e used is that measured for zero attractions, 
showing that the shear thickening stress term is independent of field strength. For panel c the value of e = is used which is 
measured at the highest packing fractions where shear thickening can be observed, showing that the phase boundary is equal 
to the shear thinning stress thb in the limit of e = 0. Solid black lines: the measured stresses at the upper and lower phase 
boundaries in the limit of zero field and small <j>. These coincide with the measured phase boundaries for B = and E = 0. 
Solid blue line: prediction of the ER yield stress from two-particle interaction (see Suppl. Mat.). Dotted blue lines: guides 
to the eye for the phase boundary between shear thinning and jammed regimes, (o): predicted r m in cases where no shear 
thickening regime was found using model predictions for j m . In each case, these values are close to or above the upper stress 
boundary, showing that the reason shear thickening was not found was because r m exceeded the shear thickening stress range. 
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FIG. 5: Fit of a stress/shear- rate curve broken up 
into shear thinning and thickening components. Data 
shown for glass spheres at E = 12.5 V/mm. Dashed line: 
fit of Herschel-Bulkley model (Eq. II D|) to the lower shear 
thinning regime. Dotted line: term ~ 7 1 ' 6 representing the 
shear thickening regime. Solid line: sum of dashed and dotted 
lines fEq. II F[) . Vertical lines: shear rates corresponding to the 
boundaries of the shear thickening range in the limit of E = 0. 



matical model, this is a good approximation as long as 
there is a sharp upturn in T(j) which is the defining fea- 
ture of 'discontinuous' shear thickening, and the shear 
thickening stress term is independent of the shear thin- 
ning mechanism. This interpretation remains true for 
e > with corrections according to Eq. II Fl It also holds 



in cases where there is a Newtonian regime before the 
onset of shear thickening, regardless of the value of e (see 
Suppl. Mat.). 

Our simple model predicts the shear thickening phase 
boundaries without knowing detailed particle properties 
or microstructure. As long as the shear thinning mech- 
anism produces a stress term that adds linearly to the 
shear thickening stress term, all sources of shear thinning 
have the same effect of hiding shear thickening, regard- 
less of the mechanisms that produces shear thickening. 
This description in terms of stress scales is not depen- 
dent on size scale and, in principle, might be applicable 
also to 'discontinuous' shear thickening in colloidal (i.e. 
Brownian) systems. In colloids, however, different mech- 
anisms for shear thickening and thinning might become 
relevant, for exam ple a shear thinning stress term due to 
Brownian motion [19( . 

The existence of an upper threshold beyond which 
shear thinning mechanisms will overwhelm shear thick- 
ening explains why in most cases attractions completely 
eliminate shear thickening Q while for some fluids with 
weak interparticle attractions it has been reported to ex- 
ist [2l|, [22| . In typical suspensions, attractions are often 
due to particle-fluid surface tension. An example is the 
common observation that cornstarch (a hydrophilic par- 
ticle) shear thickens in water but not in hydrophobic liq- 
uids [23| . One can then ask if all suspensions will shear 
thicken once the shear thinning stresses are small. In the 
experiments reported here on a variety of suspensions 
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consisting of particles including cornstarch, glass, and 
PEG, in a variety of fluids with different density match- 
ing, modified surface properties, roughness, shapes, and 
measuring conditions, we always observed 'discontinuous' 
shear thickening. Including a variety of other suspensions 
we studied, we found no examples where the shear thin- 
ning stress was small (less than the order of 5 Pa for 
particles on the order of 10-100 /im) that did not shear 
thicken at near-sedimentation packing fractions. Induc- 
tively this suggests the phenomenon of 'discontinuous' 
shear thickening is general to all hard particle suspen- 
sions at near-sedimentation packing fractions provided 
that the shear thinning stresses are below a threshold 

B|. 

The combination of ER or MR effects with shear 
thickening opens up possibilities for the design of ficld- 
rcsponsive shear thickening fluids in dampers or impact 
absorbers [f| 0, H(| . We note that earlier suggestions pre- 
sumed that the applied fields would control the critical 
shear rate [26|, but this is only true for weaker shear 
thickening (e > 0), where both the critical stress and 
shear rate vary with field. In the limit where shear thick- 
ening becomes discontinuous (e = 0), our findings show 
that the critical shear rate is controlled by the parti- 
cle packing fraction [T3], while the critical stress can be 
tuned cither passively with particle-fluid chemistry or ac- 
tively with fields. 
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I. SUPPLEMENTARY INFORMATION 
A. Additional experimental details 

In the rheometer the torque T on the tool and its ro- 
tation lo rate were measured and converted to a shear 
stress t and a shear rate 7. Both parallel plate and Cou- 
ette geometries were used (for the parallel plate geometry 
t = 2T/(irR 3 ) and 7 = Rui/d, where R is the plate ra- 
dius and d the gap size). The shear stress and shear 
rate describe the mechanical response in a geometry- 
independent form, but we do not imply or require a linear 
flow profile. Data for Figs. 1 and 3 were taken with in- 
creasing controlled stress to resolve the steep shear thick- 
ening, while data for Fig. 2 were taken with controlled 
shear rate to allow for a better fit of the Herschel-Bulklcy 
model, Eq. 1. Care was taken that no fluid extended 
outside the parallel plates and the particles were con- 
fined to the space between the plates by surface tension. 
Samples were pre-sheared for 200 s to stresses above the 
shear thickening region immediately before experiments 
commenced after which measurements were found to be 
reproducible within a typical variation of about 10-20%. 
Measurements reported were mostly taken at ramp rates 
of 500 s per decade of controlled stress or shear rate. In- 
creasing as well as decreasing ramps with different ramp 
rates were used to check for hysteresis, thixotropy, and 
transients. An example comparison of ramp rate and di- 
rection dependence is shown in Fig. [5] Some ramp-rate 
independent hysteresis was observed between increasing 
and decreasing stress measurements. While the magni- 
tude varied from suspension to suspension and typically 
about 20% of the viscosity in the shear thickening regime, 
the curves were never qualitatively different. 

For clarity, the data shown are for one ramp direc- 
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FIG. 6: Viscosity curves showing an example of hysteresis and 
ramp rate dependence for 100 fim glass spheres in mineral oil 
at <f> = 0.552. Solid line: 100 measurement points per decade 
of stress, 10 s per point. Dotted line: 50 measurement points 
per decade of stress, 10 s per point. Upper curves correspond 
to a decreasing stress ramp. Lower curves correspond to an 
increasing stress ramp. 

tion. We checked for reversibility by shearing suspen- 
sions in the shear thickening regime and then immedi- 
ately ceasing shear; the result was that the stress re- 
laxes to the zero-shear limit within seconds. Different 
gap sizes between 0.5-1 mm were used to check for fi- 
nite size effects. Reported experiments were done with 
smooth plates. Rough plates were also used in some ex- 
periments to check for slip and no significant differences 
were found. Reported packing fractions <f> are based on 
measured particle and fluid quantities mixed together be- 
fore shearing. 

To observe the discontinuous viscosity curves as e — > 0, 
confinement of the sample is important. This can be 
achieved by either using a Couette geometry or avoid- 
ing slop in a parallel plate geometry [l|, 0] . Non-density 
matched samples were measured in a parallel plate geom- 
etry to minimize the weight on the packing which pro- 
duces a yield stress in a Couette geometry [3| ■ Attempts 
to measure shear thickening of glass spheres in mineral oil 
in a Couette geometry resulted in a large yield stress due 
to this sedimentation and no shear thickening. This ob- 
servation is perfectly consistent with our conclusion that 
a yield stress from any source can hide shear thickening. 

The glass spheres were obtained from MoSci corpo- 
ration (Class IV). They were sieved through mesh sizes 
-120+170 and were measured to have a mean diameter of 
89 /xm with a standard deviation of 12 /im. For the sur- 
face tension experiments, the glass spheres were mixed 
into water after the surfactant so the total fluid volume 
matched that of the case without surfactant. This en- 
sured that the surfactant diffused throughout the sample 
and the volume fraction did not vary between the two 
experiments. These measurements were done with a par- 
allel plate setup with a 50 mm diameter rotating top 
plate with a 0.83 mm gap. 

Cornstarch was chosen as a prototypical shear thick- 
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ener for the packing fraction dependent experiments. 
Argo cornstarch was used at ambient conditions of 23°C 
and 42% humidity which included some water weight. 
The suspensions were density matched for the reported 
experiments so that the yield stress in this case was due 
to confinement Q. For the cornstarch data shown we 
used a Couette geometry consisting of a 26.6 mm diam- 
eter cylinder in a cup with a gap of 1.13 mm. We found 
that mismatching the density for starch did not affect 
the critical packing fraction because in this case the hy- 
drostatic pressure from the weight of the packing is still 
much lower than the yield stress above the critical pack- 
ing fraction. The same behavior was also obtained in a 
parallel plate geometry. 

For electrorheology measurements, any dielectric par- 
ticle in a non-conducting fluid will work. We used 
hydrophobically-coated glass in mineral oil so the 
particle-fluid surface tension was minimized. To apply 
the dc electric field, 50 mm diameter parallel metal plates 
with a gap of 0.88 mm were used as electrodes. The re- 
ported electric field value is the applied voltage divided 
by the gap size. The rotating upper plate fixture made 
electrical contact via a wire brush. This added a constant 
friction corresponding to about 0.1 Pa which limited the 
stress resolution of those measurements. After subtract- 
ing this offset, the data in the limit of vanishing applied 
field matched the zero-field value measured without us- 
ing the wire brush. Therefore, the stress resolution limit 
did not artificially set the measured onset of shear thick- 
ening. 

Most magnetorheological fluids have a yield stress even 
in the absence of a field. To obtain a sample that showed 
both shear thickening and a magnetorheological response 
we engineered a suspension to minimize particle-fluid sur- 
face tension with particles that could be filled with mag- 
netic material. To achieve this we used the PRINT® 
process. Typically, the monomer solution was prepared 
as follows: 0.30 g of magnetite (black iron oxide, av- 
erage particle diameter = 0.2 (im, Polysciences, Inc.), 
0.02 g of 1-hydroxycyclohexyl phenyl ketone (HCPK, 
Aldrich), and 0.01 g of fluorescein o-acrylate (Aldrich) 
were placed into an Eppcndorf tube followed by the 
addition of 0.1 ml of N,N-dimcthylformamidc (DMF, 
Aldrich). The monomer mixture was then mixed thor- 
oughly by vortex mixing to dissolve the HCPK photoini- 
tiator and the fluorescein o-acrylate fluorophore. Lastly, 
0.67 g of ethoxylated(20) trimethylolpropane triacrylate 
(MW = 1176 g/mol, SR415, Sartomer) was added to the 
monomer mixture and vortex mixed again. The result- 
ing solution was composed of 30% (w/w) magnetite, 67% 
(w/w) triacrylate, 2% (w/w) HCPK, and 1% (w/w) flu- 
orescein o-acrylate. The rod-shaped particles were then 
fabricated using the PRINT process, which has been de- 
scribed elsewhere [U, 0] ■ Molds for fabrication of PRINT 
particles were supplied by Liquidia Technologies. For 
the magnetorheological experiments, the particles were 
suspended in poly(cthylenc glycol) dimethyl ether (Mn 
= 500 g/mol, Aldrich). These measurements were con- 



ducted in a 20 mm diameter parallel plate geometry with 
a gap of 0.9 mm. For imaging purposes, DyLight 549 
Maleimidc (MW = 1007 g/mol, Fisher) was used as the 
fluorophore. 



B. Yield stress 

The yield stress can be defined differently and thus 
measured in several different ways. A static or dynamic 
yield stress can be measured for either increasing or de- 
creasing control ramps, respectively, and each can be 
done with either controlled stress or shear rate. In our 
experiments, each method yielded similar yield stress val- 
ues. Some hysteresis was observed between the static and 
dynamic yield stresses, which was larger for faster ramp 
rates. At slower ramp rates the hysteresis loops con- 
verged to a relatively small difference (less than a factor 
of 2). The reported data were taken at ramp rates in this 
latter regime. 

By defining the viscosity as rj = r/7, it is infinite be- 
low the yield stress since the shear rate is zero. Given 
that shear thickening requires the viscosity to increase 
with stress, ?/ must first drop to finite values, so a viscos- 
ity function with a continuous first derivative necessarily 
displays shear thinning before entering a shear thicken- 
ing region. A different value for the yield stress does not 
change this general behavior, but can move the onset of 
the shear thickening region. The shear thickening stress 
range can therefore depend somewhat on the yield stress 
definition, specifically if the lower shear thinning region is 
small. For the purposes of comparing shear thinning and 
shear thickening stresses to determine the shear thick- 
ening regime, the particular criterion for evaluating the 
yield stress is irrelevant as long as it is done consistently. 
The conclusion that the shear thickening phase bound- 
ary is determined by the shear thinning stress does not 
depend on which yield stress is measured or any specific 
form for the model. We chose the Herschel-Bulkley model 
with exponent 1/2 only because it fits the data well (but 
see below for other exponents). 



C. Connection between particle interactions and 
macroscale rheology 

The connection between field-induced interparticle at- 
tractive forces and the yield stress can be explained 
through electrorheology models Q . In an applied electric 
field E the induced dipolc moment density scales as [3eoE 
where f3 is an effective dielectric constant that depends 
on particle and fluid dielectric constants and saturates 
at values of order unity for all but a near-exact dielec- 
tric match. The resulting net force between neighboring 
particles scales as F ~ eo^LP 2 E 2 a 2 for particle diame- 
ter a and liquid dielectric constant e^. The yield stress 
can be obtained by dividing this attractive force by an 
effective particle surface area, giving a yield stress scale 
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T y = 12ne e L (3 2 E 2 Q. This result is shown in Fig. 2b and 
agrees with the measured yield stress at high field values. 
Settling becomes more important below the gravitational 
stress scale Ap<?a w 1.5 Pa where Ap is the density differ- 
ence and g is the acceleration of gravity. If the attractive 
stress is below this scale, particles will tend to settle at 
lower shear stresses rather than form chains to jam the 
system. This is likely the reason that the yield stress falls 
below the E 2 scaling at lower field values. The agree- 
ment of the yield stress with the attractive stress scale 
shows that the yield stress and hence the shear thicken- 
ing regime can be connected to individual particle prop- 
erties. It is interesting to note that this calculation did 
not require any specific knowledge of the flow or packing 
structure. 

The yield stress scale could be put in terms of the at- 
tractive force divided by particle surface area, or cquiv- 
alcntly the attractive energy per unit volume. This at- 
tractive stress scale can be calculated for other types of 
attractions as well to relate the yield stress to microscopic 
properties. For example, we can estimate the expected 
yield stress from other sources of interactions that might 
be operative between cornstarch particles. To check for 
this, we used optical tweezers to place two cornstarch 
particles next to each other in water and allowed them 
to diffuse. An attractive or repulsive potential can be 
measured by observing the probability distribution of the 
separation distance over time. In the tweezer experiment 
the resolution was about 1 pN and down to this instru- 
mental limit no attractive or repulsive forces were ob- 
served. Dividing this value by particle surface area puts 
the upper limit on the yield stress due to attractions at 
around 10 -2 Pa. This is consistent with the fact that 
we did not observe any yield stress in cornstarch suspen- 
sions down to our instrument resolution of 10~ 3 Pa at 
low packing fractions. 



D. Approximations of the phase boundary 

To understand how the lower shear thickening phase 
boundary is determined by the shear thinning stress, 
we next discuss various approximate solutions based on 
Eq. 3. The basic idea is as follows: Given that the shear 
thickening stress is independent of the strength of at- 
tractions (see Fig. 2) and described by Eq. 2, the phase 
boundary can be determined by measuring j m , a 2 and e 
for zero field, and then calculating t v (B) and a\{B) as a 
function of applied field. To check the feasibility of this 
approach, we first show that the fit parameters ai and 
e are independent of the applied field. This is done by 
fitting the experimental data to Eq. 2 for various field 
values, covering a range up to 80% of the shear rate at 
the viscosity maximum or up to 3 Hz if there is no max- 
imum (this covers roughly the same shear rate range). 
The result, shown in Fig. confirms the claim that the 
shear thickening stress is unaffected by (field-induced) 
attractions. For brevity, we will show the subsequent 
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FIG. 7: Parameter values obtained from fitting Eq. 2 to the 
data on shear thickening with the MR effect, (o): 0,2 in units 
corresponding to Pa for stress and Hz for shear rate. (•): 
1/e. No resolvable trend in either parameter is found as B is 
varied. 



analysis only for the magnetorhcology data but it applies 
for the electrorheology data as well (the <j> dependence is 
more complicated because the packing fraction controls 
the shear thickening behavior itself in addition to tuning 
the yield stress). 

The simplest approximation of the phase boundary is 
to assume the discontinuous limit e = 0. When varying 
the yield stress via attractive interactions, the fact that 
the shear thickening part of r(j) is steep means the the 
onset of shear thickening occurs at a nearly constant 7 m . 
This allows for a simplification since we can approximate 
7 m in Eq. 3 by its value measured for zero applied field. 
This gives r m = THBijm) indicating that the stress at 
the phase boundary is equal to the shear thinning stress. 
We evaluate Eq. 1 using fit values of r y (B) and ai(B) for 
each applied field evaluated at 7 m ,0j the measured onset 
at zero applied field. This is shown as purple symbols in 
Fig. [5J along with the data from Fig. 4a. This e = ap- 
proximation underestimates the onset of shear thickening 
by 50 ± 20% (errors indicate a standard deviation). The 
fact that this approximation gives the threshold where 
attractions begin to move the onset and the increase in 
the onset with field within about a factor of 2 confirms 
that the phase boundary is determined by the shear thin- 
ning stress. 

A better quantitative match to the lower phase bound- 
ary can be obtained by accounting for the non-zero e. 
The orange symbols in Fig. [5] are plotted for the same 
T y and a\ as before but now using the measured value 
of e = 0.55. This increases the predicted r m by a factor 
of 1.5. As a result, the measured boundary is underes- 
timated by 9 ± 9% in the low-field region and at higher 
fields, where attractions are reducing the shear thicken- 
ing regime, by 33 ±20%. Overall, this better predicts the 
point where the attractions are strong enough to increase 
the onset of shear thickening but still underestimates the 
effect of attractions. 

The next correction is to account for the change in j m 
with attractions for e > 0. Using techniques similar to 
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FIG. 9: Solid triangles (T): measured shear rate at the onset 
of shear thickening 7 m . Open cirlces (o): calculated jm from 
Eq.g] 



FIG. 8: Phase diagram showing different approximations of 
the shear thickening phase boundary. Lower (T) and up- 
per (A) boundaries of the shear thickening regime are as in 
Fig. 4a. Purple circles (•): estimate in the limit of e = 
giving Tm = THB^m) evaluated at j m ,o. Orange circles (•): 
estimate accounting for the non-zero e using the measured 
value of e = 0.55 in Eq. 3. Open red circles (o): estimate 
further accounting for the change in j m with e and B by eval- 
uating Eq. 3 at j m calculated from Eq. [4] Solid red circles 
(•): Eq. 3 at the measured y m . Solid green circles (•): pre- 
diction from Eq. 3 at 7 m calculated from Eq. [4] using only 
data obtained at B — and the lower shear thinning regime 
for larger B. Open green circles (o): evaluating Eq. 3 at j m 
calculated from Eq.[8]with a = 1/2. Dotted green line: fit of 
the open green circles indicating the phase boundary between 
shear thinning and Newtonian regimes. 



those used in Sec. II Fl an exact implicit equation can be 
written for j m in a form that shows how j m varies with 
r: 



im{BY 



(7m,o)" 



2(1 -e)aa 



Aai(B) 



2r y (B) 



where Aai(B) = a\{B) - a±(0). 



V^iB) 
(4)' 

Eq. H] reduces to 



7m = 7m, o for no attractions [Aai(B) = 0, r y {B) = 0] 
as expected or for step-function stress-shear rate curves 
(e = 0) as already claimed, and it justifies the simpli- 
fication T m = THB(jm,o) in the limit of e = 0. Such 
constant onset shear rate when the yield stress is varied 
in the limit of e = contrasts with the constant onset 
stress when the packing fraction is varied (not including 
the contribution of the yield stress) P, 0, @ • 

Since Eq.|4]is an implicit equation, it must be evaluated 
numerically. We note that since the yield stress is a small 
contribution to the overall shear thinning stress as seen in 

1/2 

Fig. 4a, i. e. ai7m 3> r y , the rightmost term with j m (B) 
in the denominator is small. Starting with the value of 
7 m (_B) = 7m, o °n the right hand side, we can evaluate 
Eq. @] iteratively. The value of j m converges to within 



a few percent after only 2 iterations. Thus for a simple 
explicit estimate one can set j m (B) = A / m ,o on the right 
side. This estimate of 7 m is shown in Fig.[9]in comparison 
to the measured j m . It is seen that the model captures 
the increase in j m i n with attractions, particularly the 
point where attractions start to increase j m which occurs 
at the same point the stress starts to increase due to 
the yield stress pushing up the onset of shear thickening. 
Beyond that point the model overestimates the measured 
values by 30 ± 20%. 

We now evaluate Eq. 3 using the calculated value of 
j m (B) from Eq. @] e = 0.55, and the fit values of t v (B) 
and a\(B). This is shown as the open red symbols in 
Fig. [5] (same as in Fig. 4a) and gives the entire phase 
boundary within a standard deviation of 28%. For the 
packing fraction dependence in Fig. 4c, there is no com- 
parable prediction for the onset shear rate because the 
shear thickening term varies with packing fraction. Thus 
the open red symbols in panel c correspond to Eq. 3 eval- 
uated at the smallest measured onset shear rate. 

A check on the validity of Eq. 3 for describing the phase 
boundary can be made by using the fit parameters t v (B) 
and a\{B) and j m (B) measured at the minimum of rj(r). 
This is shown as solid red symbols in Fig. [8] (same as in 
Fig. 4a) which agrees with the measured phase boundary 
to within a standard deviation of 12%. For comparison, 
if we repeat measurements keeping all control parameters 
constant, the typical variation in the measured r m is 11%. 
Thus, the model is accurate in describing the onset of 
shear thickening up to the resolution of the data. 

Given the assumptions that the shear thinning and 
thickening terms add linearly and the shear thickening 
term is independent of attractions, in principle we can 
predict the phase boundary using only the shear thick- 
ening curve at zero field and the effect of attractions on 
t~hb- The above analysis was all done using fits of data 
up into the shear thickening regime. To show the pre- 
dictive power of the model, we now perform the fit to 
the data for zero field only to obtain 02, e, and 7 m ,o an d 
keep these fixed. We then fit Eq. 2 to data for non-zero 
field only up to some cut-off 7 < 0.3 which is in the lower 
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shear thinning region for data with non-zero field. This 
allows us to obtain r y and ai . The fitting cutoff can be 
chosen based on the zero-field data because the attrac- 
tions always increase the onset of shear thickening. We 
then evaluate Eq. 3 using the fit values of r y and eti, and 
7 m from Eq. 21 This gives the solid green symbols shown 
in Fig. [Sk and b. 

This prediction of the phase boundary without using 
any data from the shear thickening regime except at zero 
field overestimates the phase boundary by 20 ± 60%. 
While this agreement is not as good as when we fit the 
full data set, it shows that the effect of attractions on 
shear thickening can be predicted within about a factor 
of two. In the case where packing fraction is varied, the 
shear thickening stress itself changes; a 2 and e vary with 
c/>, which shifts j m without the influence of attractions. 
Thus cither a 2 and e, or j m must be obtained as a func- 
tion of 4> to determine the phase boundary. 

A Newtonian regime is sometimes found before the on- 
set of shear thickening, for example in Fig. 2a. We did 
not explicitly include a Newtonian term in the model. 
However, the generalized expression for r m in Eq. [7] can 
apply for a Newtonian stress term when a = 1 (see the 
next Section, below) . The onset of shear thickening can 
still be expressed by r TO in Eq. [7] if a Newtonian stress 
term is added to Eq. 2, although the value of 7 m ,o would 
generally increase. We expressed Eq. 0] in a perturbativc 
form rather than as a simple dependence on the fit pa- 
rameters so it still applies in the case where there is a 
Newtonian regime whether or not it can be described by 
a linear addition to the stress-shear rate relation of Eq. 2. 
If instead Eq. [5] is used for j m without consideration of 
a Newtonian regime, the phase boundary would be un- 
derestimated for weak attractions (open green symbols 
in Fig. Thus the difference between the open red and 
green circles in Fig. [5] is due to the Newtonian regime. 

For simplicity we omitted any Newtonian regime from 
the main paper and included it in the shear thinning 
regime. This does not change the conclusions but to be 
more general we can restate them in a way that includes 
the possibility of a Newtonian regime. This regime dis- 
appears for stronger attractions when the shear thinning 
stress overwhelms the Newtonian stress term. Thus for 
shear thickening to occur in general, the shear thicken- 
ing stress must overcome the sum of shear thinning and 
Newtonian stresses. On the other hand, for attractions to 
affect the onset of shear thickening, they must exceed a 
threshold equal to the inherent shear thinning and New- 
tonian stresses at the onset. 



E. Notes on possible mechanisms for shear 
thickening 



have successfully described Continuous Shear Thickening 
which occurs at lower packing fractions and higher shear 
rates, they have not been able to reproduce the steep 
stress/shear-rate relation (e « 0) characteristic of Dis- 
continuous Shear Thickening, instead the smallest value 
of e allowed in those models is 1/2 M. Inertial granu- 
lar models have a similar limitation [lfjj. In addition, 
hydrodynamic models predict that the shear thickening 
stress should be affected by attractions [ll|, E3- How- 
ever, for Discontinuous shear thickening this is not what 
we observe. 



F. Derivation of Eq. 3 

Here we derive an expression for the onset of shear 
thickening from t(7) given in Eq. 2 with a generalized 
Herschel-Bulkley form for the shear thinning term thb = 
T y + ai7 Q . The onset corresponds to the local viscosity 
minimum which satisfies 



cry 



Rearrangement gives 



(5) 



a 2ill" = + t 1 ~ a ) 0l ^n] ■ ( 6 ) 

Substituting Eq. |6] into Eq. 2 evaluated at j m gives 



T m = T HB (7m) + \- T V + a ) a l7m] ■ (7) 

Setting a = 1/2 gives Eq. 3. This shows that in the limit 
of e = the onset of shear thickening is equal to thb (7m) 
regardless of the form of the shear thinning term. 



G. Derivation of Eq. [4] 

Here we derive the expression for the shear rate at the 
onset of shear thickening 7, n , similar to the derivation for 
T m (7 m ). We rearrange Eq. [5] to get 



, - —a 

7m 



(l-e)oa 



[(1 - a)a x + rj,7 ro 



(8) 



To put this in the form of Eq. [4] to directly describe 
the perturbation in j m with an additional shear thinning 
term, we evaluate this for zero additional attractions to 
obtain 



The phenomenological approach presented here does 
not address the microscopic origin of shear thickening in 
suspensions, but the data put constraints on the region of 
validity for existing models. While hydrodynamic models 



7m, 



(1 - ()a 2 



{I -a) (ai-Aai) 



(9) 



where Aa± = a± — aifi and ai.o is the value of a\ for the 
unperturbed state without the additional shear thinning 
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term. Substituting this back in to Eq. [8] to eliminate a\ 
results in 

W° = + ( , 6 s [(1 - a)Aai + T yl - a ] . (10) 
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